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THE RESPONSE OF A SINGLE WALL SPACE STRUCTURE TO IMPACT 
BY COMETAEY METEOROIDS OF VARIOUS SHAPES 

By 

Robert J. Hayduk 
ABSTRACT 

Linear, small- deflect ion plate theory is used to study the stress 
at the contact axis and the deflection of an infinite plate caused by 
the impact of an axisymmetric coraetary meteoroid. The analysis assumes 
that momentum exchange is the primary mechanism, that the time of exchange 
is instantaneous, and that the momentum of the meteoroid is negligible 
after impact. The stress at the origin is reduced to a single definite 

: *r. ; 

integral and the deflection to the Hankel inversion Integral, both 
requiring definition of the particular projectile before further evalua- 
tion. A particular cometary meteoroid is mathematically represented in 
the analysis by its projected momentum per unit area onto the plate. 

The three specific shapes studied are the usual projectile shapes 

* 

used in hypervelocity laboratories - cylinder, cone, and sphere - even 
though the analysis is not intended for the high-strength, high-density 
laboratory projectiles. Projectile comparisons based on equal mass, 
diameter, and total momentum indicate that frangible, low-strength cone 
projectiles cause significantly higher stresses and larger displacements 
of the plate at short times after impact than s imilar sphere and cylin- 
der projectiles. 
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Y. INTRODUCTION 


The literature indicates general concern over the hazard of travel 
in the meteoroid environment of space. Catastrophic failures of manned 
and unmanned space vehicles can result from meteoroid penetration of an 
underdesigned spacecraft hull. The National Aeronautics and Space 
Administration is investigating the reactions of space structures sub- 
jected to impact by hypervelocity projectiles. These investigations 
along with added knowledge of the meteoroid environment are continually 
upgrading the spacecraft design equations. 

Present laboratory simulation capabilities are barely reaching the 
lower limit of the meteoroid velocity range with moderately large projec- 
tiles. Since space structures cannot be subjected to an extensive simu- 
lation of the environment , any design equation can only receive limited 
verification before being extended to the environment. Any analytical 
model, then, should be as realistic as possible in the approximation of 
the hypervelocity phenomenon of meteoroid impact. 

This thesis presents an analytical model of the hypervelocity impact 
between an axi symmetric cometary meteoroid and a single wall space struc- 
ture. Linear, small -deflect ion plate theory is used to study the stress 
at the contact axis and the deflection of the single wall. 

Classical plate theory has been used by many investigators to study 
the response of plates to projectile impact. Reference 1 contains solu- 
tions to impulsively loaded infinite plate problems frequently referenced 
in the literature. Particular problems are an impulsive point force at 
the origin, an impulsive force uniformly distributed over a circle, and 
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an impulsive Gaussian distribution of pressure. In reference 2, the 
deflections away from the impact point of a large plate centrally 
impacted by a bullet agreed after a short time with those obtained by 
mathematically simulating the Impact by an impulsive point force of 
finite duration on an infinite plate. Also presented in reference 2 are 
solutions to the impulsive point force problem for various time func- 
tions: finite duration, step function, and unit impulse function. 

The hyper velocity impact between a right circular cylinder and a 
thin elastic plate was analyzed in reference J. It was assumed that the 
projectile gave the contact portion of the plate a uniform initial 
velocity determined from conservation of momentum. The stress at the 
axis of symmetry was used as an indicator of the failure threshold to 
obtain a simple ballistic limit thickness expression. The analysis pre- 
sented here was performed using an approach similar to that used in 
reference 3, but generalizes the problem to any axisymmetric projectile. 
The displacement of the single wall and the time variation of the stress 
at the origin are investigated. Solutions for three particular projec- 
tiles - cylinder, cone, and sphere - are presented for comparison. 


VI. SYMBOLS 


a 

B(x,y) 

c 

D 

E 

e 

P 

g 

h 

J 


radius of contact betveen projectile and plate 
beta function 

speed of sound in plate material, VE/pt 

flexural rigidity of the plate, Eh3/l2(l - v 2 ) 

elastic modulus of plate material in tension and compression 

strain 

hypergeanetric function 
projectile functional 
target thickness 
Bessel function 


K a 


2a 2 

ch 


- v 2 ) 


M 

Z 

P 


bending moment resultant 
projectile length 
Hankel transform parameter 


Pt) v b 

« > K 

p t h 


r 

t 


n 

V 

p 


r 

a 


radial coordinate 
time 


V velocity 

w(rj,t) transverse displacement of the plate's midplane 


Poisson's ratio 
density 
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a stress 



Subscripts 
t target 

p projectile 

co cone 

cyl cylinder 

r radial 

sph sphere 

9 tangential 

p indicates order of Bessel function 

v indicates order of Bessel function 

Superscripts 

• derivative with respect to time 

~ nondimens lonal 

- transform 

» nondimensional equl-raass and -diameter 

|i indicates order of Bessel function 

v indicates order of Bessel function 



VII . MATHEMATICAL DEVELOPMENT OF THE PROBLEM 


Transformed Governing Equation and General Solution 
The governing equation for the bending of a thin plate from linear, 
elastic, small-deflection theory is 

V v » -K? gjjf (la) 


where 


K - ir - *“> <»> 

and the operator in polar coordinates is 



w( T] ,t ) is the transverse displacement of the midplane of the plate, 
a is the impacting projectile radius, c is the speed of sound in the 
plate material, h is the plate thickness, v is Poisson , s ratio, and 
t| is the radial coordinate r divided by the projectile radius . 

Having conveniently selected the plate as infinite in extent makes 
the problem amenable to a zeroth order Hankel transform in the nondimen- 
sional radial coordinate T). The forward and inverse transforms are, 
respectively, 

n oo 

f(p) = / T)f(n)J 0 (pt])dT) (2a) 

J o 
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(2b) 


f(n) s / p^Cp)*©^ 1 ^ 

w 0 


Multiplying the governing equation (la) by r>X 0 (pn) and integrating with 
respect to rj by parts several times yields the following result: 


:v 0 <w> r. 


1 <3 / c3w' 

n §n, 


+ p 


a / <w 

3n v a 


^(Pl) 


w 


- p 2 n J 0 (pn) - p^w^Cpn) 


+ kQ ^5 ^ Ji® 


\ 


J 


dt 2 


+ pV a Q 


(5) 


T)=0 


If the boundary conditions of zero deflection, slope, and shear per 
unit length at infinity and zero slope and finite shear at the origin 
are assumed and coupled .fith the kernel of the transformation, the term 
in brackets vanishes. If the time is nondimenslonalized to r » g, the 
transformed differential equation becomes 


d 2 w 

dr 2 


+ p^w = 0 


with general solution 


(4) 


w(p,r) = A(p)sin(p 2 T ) + B(p)cos(p 2 T ) (5) 


Applying the zero deflection initial condition (w(p,0) « 0) indicates 
that B(p) must be zero. The coefficient A(p) is determinable from the 
transform of the initial target velocity, that is, 
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Therefore, 



1 dw(p,T) 


sip 2 A(p)cos(p 2 t) 


A(p) a Kp” 2 w(p,0) 


( 6 ) 


Determination of Initial Plate Velocity 
(Cometary Meteoroid Impact Model) 

Reference 4- summarizes the interplanetary meteoroid environment as 
consisting of debris in two distinct forms: hard, dense asteroidal 

particles in relatively sparse quantities and porous, fragile cometary 
particles constituting almost 100 percent of the total. These cometary 
meteoroids, as described in reference 5, are thought to have a density 
on the order of 0.5 gm/cm^ and extremely low crushing strength such as 
cigar ash. 

Upon impact with a structure at a velocity between 11 and J2 km/sec, 
a cometary meteoroid is assumed to act as a projectile consisting of 
discrete small particles which interact negligibly since the crushing 
strength is so low. This leads to a model of continuous momentum 
exchange during impact if the further assumption is made that a momentum 
balance is of primary importance. This time variation of the loading, 
so to speak, is extremely difficult to include in an analysis. There- 
fore, this model will assume the worst case - all the momentum of the 
cometary meteoroid impacts Instantaneously. 
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Any axisymmetric projectile will be represented in the analysis by 
the function representing its total momentum projected onto the plate 
rather than the analytic expression for its shape. Figure 1 illustrates 
this point for a sphere. The upper diagram shows a sphere and its 
analytic expression 


z 




The lower diagram shows the functional representing the total projected 
momentum per unit area of the sphere, that is, 


^sph^) “ 2a 



( 8 ) 


Hence, the momentum Impacting an annular ring of radius T) and 
width dt) is 


2*a 2 n if) f p (T))Vpp p (9) 

Balancing momentum between the projectile and target yields the initial 
velocity of the target as 


w(t],0) = 


V p Vv \ 

\~K m tr) f -n( r \)f o < t) < l 

\P t Ah/ P 


1^0 ; 1 < tj < » 


( 10 ) 


where the momentum of the projectile has been neglected after impact. 
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z 



Momentum Per Unit Area Functional: fspt/ 7 ?) ■ 2a(l - 17 2 ; 


Figure 1.- Comparison of the sphere's geometric equation and 
functional form of the projected momentum per unit area. 
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Transforming this initial velocity evaluates the coefficient A(p) 
(eq. (6)) as 


a(p) = Q z pP ~ 2 f ngp(n)J 0 (pn)dri 

= Q^pP" 2 gp(p) 


( 11 ) 


P v 

r-n T 


vhere the nond^ns tonal entity K ^ has been represented by Q, 
and the functional fp( tj) has been represented by a characteristic 
length, Zp, times another functional gp(il). 



VIXI . MATHEMATICAL SOLUTION OF THE PROBLEM 


Displacement of the Plate 

Substituting equation (ll) and B(p) » 0 Into equation (5) and 
inverting according to equation (2b) yields the displacement of the plate 
in integral form as 


n oo 

Wp( n^T) ^ Q^ p J P^^pJJqCpt) )sin(p 2 T)dp (12a) 


or dividing by QZp to nondimensionalize yields 


Wp(rjrt) = f P^gpCpJJoCp 1 )sln(p 2 ' r )dp 

v r\ 


(12b) 


Stresses, Strains, and Bending Moment Resultants 

at the Origin 

From linear, elastic plate theory, the appropriate equations for 
stress, strain, bending moment resultants, and shear stress resultant are, 
respectively, 

o. 


'rr 


zE /c^w v 
a 2 (l - 7)W + 


(15a) 


zE 


7 00 = 


1 + v d 2 w 

a 2 (l - v 2 ) ^ dt] 2 


-rr 


z d 2 ^ 

2? Sip 


( 15 b) 


(14a) 


z 1 

e e0 = - -s 3rj 


(i4b) 


li 
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Hj. 


Mq « - 


D ) bfiw . V 






+ -r 



1 chf , 

n Stj 


afw 


(15a) 

(15b) 



- £ SJ^ 


(16) 


where the remaining quantities are zero because of thin plate assumptions 
and the axisymmetry of the problem. 

Due to symmetry and the monotonlcally decreasing projectile shape the 
maximum stress, strain, and bending moment resultants occur at the axis 
of symmetry. Knowledge of these maximum plate responses is essential 
for hypervelocity impact studies of spacecraft with protection of the 
occupants (both equipment and people) in mind. 

At the axis of symmetry (origin of coordinate system), both the 
radial and circumferential stresses, strains, and bending moment resul- 
tants are identical because 


11m 
T) — > 0 


1 bw d^w 

“ SJf 


T)=0 


( 17 ) 


by 1 'Hospital's rule. Thus, the nondlmensional stress, strain, and 
bending moment resultants are identical at the origin, that is, 


a 


/v# 

e = M = 




n=o 


(18) 
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where 


y a a2 .(i zjl 

zEQZp 


cr rr or 0Q0 


TJ»0 


(19) 


~ a' 
e a 


z QZ 


P 


e rr or e 00 


n*o 


M 


a‘ 


D(1 + v)Q2 p 


M r or M0 


T )=0 


( 20 ) 


( 21 ) 


Before differentiating the displacement (eq. (12b)) to evaluate the 
stress at the origin, substitute the integral expression (eq. (ll)) for 
gp(p) and switch the order of Integration to get 

n 1 n oo 

Wp(n,T) = J ugp(u)du J p" 1 J 0 (pu)J 0 (pT))sln(p 2 T)dp (22) 

Differentiating twice with respect to T) and letting T) go to zero 
gives 


n2~ 

o w- 


P 


8t) 2 


TJ=0 


= - | J ugp(u)du J pJ Q (pu)sln(p 2 T )dp (25) 


and the stress at the origin (eq. (l8)) can be written as 


dp = “ J ugp(u)du J pJ 0 (pu)sin(p2r )dp 


(24) 
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After evaluation of the Infinite Integral (see appendix), the stress 
at the origin reduces to a single definite integral, that is, 


®P “ 57 / o u«p(u)c0B^au; T / 0 


(25) 


which can be integrated once gp(u) is known for the particular axisym- 
metrlc projectile. 
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IX. PARTICULAR SOLUTIONS FOR CYLINDRICAL. CONICAL. 


AND SPHERICAL PROJECTILES 

The three particular axlsymmetrlc projectiles to be studied are the 
usual ones used in hypervelocity laboratories - cylinder, cone, and 
sphere - even though the analysis is not Intended for the high-strength, 
high-density laboratory projectiles. Since it is assumed that the total 
momentum of the projectile is imparted to the plate instantaneously, 
this analysis does not distinguish between the conical projectile (or 
any other projectile for that matter) impacting base down from the same 
projectile impacting apex down. 

Figure 2 presents the three particular axlsymraetric projectiles and 
their functional representation in the analysis. 

Displacement of the Plate 
Cylinder . - g cyl (n) « 1 

Substituting 6 C yi(il) - 1 into equation (ll) yields 

nl 

g C yl(p) s J uJ 0 (P u ) du 

» p" 1 J 1 (p) (26) 

The nondimens lonal displacement (eq. (12)) for the cylindrical projectile 
thus becomes 

/ oo 

P~ 2 J!(p)J 0 (pT])sln(p 2 T)dp (27) 

o 


15 
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™_2a— • «. 

2 a * 


Projectile : 

Cylinder 

Cone 

Sphere 

1 

f p (7 7>“t 

l c y|( 1 ) 


dm 

2a(l-7? 2 ) 2 


Figure 2.- Specific projectile shapes studied and their mathematical 

representation in the analysis. 
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cone.- g co (n) ® 1 - n 

Substituting g co (r|) 3 1 - H Into equation (li) yields 



(p) uJ Q (pu)du - 


u 2 J Q (pu)du 


* P’ 1 Ji(p) 


1 v ( l . .1 
3 r2\2' X '2' 



(28) 


(see appendix for evaluation of second integral) which when substituted 
into equation (12) gives the nondimensional displacement of the plate 
for the conical projectile as 


w C 0 (TbT) = f p“ 2 J 1 (p)J 0 (pTl)sin(p2r)dp 
J o 

- | f P" 1 f P 2 ) J o^ pTl ^ sin ^ p2T ^ dp 

(29) 


The first integral represents the deflection of the plate due to a 
cylindrical projectile of the same length as the cone, and the second 
integral represents the correction necessary to account for the cone’s 
taper. 

Sphere .- g sph (T)) 3 (l - T ) 2 ) 1 / 2 
Substituting g sph ( tj ) = (l - T) 2 ) 1 ^ 2 into equation (ll) yields 




J Q (pu)du 



( 50 ) 
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(see appendix for evaluation of this integral) which when substituted 
into equation (12) gives the nondimens ional displacement of the plate 
for the spherical projectile as 

^sph^' 7 ) * (|) ^ f P" 5/ ^ J3/ 2 (p) J o(P T l) siri (P 2T ) d P (31) 

The sine function in the integrands of equations (27), (29), and 
(51) can be approximated for r « 3, by the first term of the sine 
series, that is, 

sin p^r « p^T 

This substitution renders the short time displacements for the three 
cases integrable (see appendix) yielding for the cylinder 


for the cone 


v cyl (n,T « 1) = 


Tj 0 < TJ < 1 
Oj 1 < T) < oo 


( 52 ) 


W co (i|,T « 1) = 


t (1 - Tj); 0 < TJ < 1 


; 1 < n < 


00 


(33) 


and for the sphere 


v sph (n,T « = 


t(1 - T \ 2 ) 1 ^ 2 i 0 < T) < 1 
0 } 1 < T) < 00 


(3*0 


19 


For very short time, then, the deflection of the plate is linear in 
time, is confined to that portion of the plate impacted, and has the 
shape of the projected momentum per unit area of the projectile. 

Tnus far, in the three examples, only projectile shape has been 
specified. To make realistic comparisons of the plate’s deflection at 
various times, select projectiles of equal mass and diameter. The 
characteristic lengths then vill be Z cyl = - a, Z co « 4a, and 
Z sph = 2a. The comparable nondimens lonal deflections can be written as 


5 n * 

<*1 2 aQ 

2 <v 

5 Vi 

(35a) 

= v co 

co " 2 aQ = 

^co 

(35b) 

a W Sph 

w sph = 2 aft 

= "sph 

(35o) 


The short time displacements of the plate for comparable projectiles 
are plotted in figure jj. Figures 4 through 7 present plate deflections 
at t = 0.15, 0.4, 0.6, and 1.0, respectively. As early as t = 0.15 
the shape of the projected momentum per unit area of the projectile is 
no longer evident in the plate's deflection. As time increases, the 
velocity of the plate origin decreases faster than at positions away from 
the origin and a flexural wave moves radially outward. 

The maximum deflection (tj = G in eqs. (52), (55), and (54) com- 
bined with eqs. (55a) , (55b), and (55c)) at a particular time (t « 1) 
due to the cone impact is three times and due to the sphere is one and 
one-half times that for the cylinder. As early as t = 0.15 these large 








2k 


0 


.5 


1.0 


w X 10 1.5 


2.0 


2.5 


3.0 



Figure 7»- Comparison of the plate's 
impacting equi-mass and -dii 
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differences In deflection have drastically diminished ( ■ » 1.57 and 

f . \ \ v cyl 

* l.l8) as shown in figure By r * 1.0 the deflection of the 
v cyl / 

plate is nearly the same for all three projectile impacts. 

The short time and r » 0.15 plate deflection curves indicate that 
the cone projectile causes the largest deflection at early times after 
impact. The sphere causes correspondingly less deflection of the plate 
than the cone, but more than the cylinder. 

Maximum Stress for Specific Cases 
Since the maximum stress, strain, and bending moment resultants 
have been shown to be of the same functional form, the discussion will 
be restricted to maximum stress with the understanding that the comments 
are applicable also to the maximum strain and maximum bending moment 
resultant . 

Cylinder . - g cyl (t)) * 1 

Recalling equation (25) and g cy i(n) = 1; the nondimens lonal stress at 
the origin for the cylindrical projectile is 


** X 

0 cy! “ ft 



(56a) 


Vi ■ I sln ( 4T) 


(36b) 


which readily integrates to 
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This result indicates that the nondimensional stress at the origin 
oscillates between ±l/2 with initially infinite frequency and diminishing 
frequency as * r increases. 

Cone . - 6 co (t)) » 1 - h 

Substituting g co (n) * 1 - h into equation (25) yields the nondimen- 
sional stress at the origin as 

■ w f* u(1 - u)oos (^) du (57a) 


which, after integration by parts and an appropriate change of variables, 
becomes 


co 


» T 1 ^ 


L 


(kr) 


-l/2 


sin y 2 dy 


( 37 b) 


The Integral is the familiar Fresnel Sine Integral which is tabulated 
in mathematics handbooks. 

Sphere .- g sph (T)) = U - rft) 1 / 2 

The sphere always presents problems in impact analyses; this one not 
being the exception. Rather than getting nice, clean expressions for 
the stress at the origin as equations ( 56 b) and (27b) for the cylinder 
and cone, respectively, the sphere gives a series answer to the integral 
resulting from substituting g sph (n) * (1 - n 2 ) 1 ^ 2 into equation (25), 
that is, 


a - A. 
sph - 57 




(38a) 
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which Is evaluated In the appendix as 



1 1,1/2 1 


(-X) n (ltT)’ (2n+:l) 



(58b) 


Again in the three examples, only projectile shape has been specified. 
Selecting the projectile characteristics as before, the comparable non- 
dimensional stresses at the origin become 


a &(l m v) ^ 2 y 

cyl “ 2zEQ cyl 81 3 a cyl 

(59a) 

C G0 * C co a 2c c0 

2zEQ 

(59b) 

=* a(l - v) ** 

a sph 2zEQ a sph ® 0 sph 

(59o) 


Equations (59a), (59b), and (59c) combined with equations (56b), (57b), 
and (58b), respectively, have been evaluated for t ranging from zero 
to 1 and are presented in figure 8 to demonstrate the effects of these 
projectile shapes on the maximum stress at the origin. Note that in all 
three cases a distinct peak occurs after which the stresses rapidly con- 
verge. The peak stress for the cone impact is 1.58 times and for the 
sphere impact is 1.22 times that of the cylinder. Hence, the peak stress 
at the origin also indicates that the cone causes larger responses of 
the plate than either the sphere or cylinder. 




x. CONCLUDING REMARKS 


Linear, elastic, small-deflection plate equations were used to study 
the response of a single wall structure to the impact of an axisymmetric 
cometary meteoroid. The analysis assumed that momentum exchange is the 
primary mechanism, that the time of exchange it instantaneous, and that 
the momentum of the meteoroid after impact is negligible. 

Three particular shapes were studied - cylinder, cone, and sphere. 
Projectile comparisons based on equal mass, diameter, and total momentum 
indicate that frangible, low-strength cone projectiles cause signifi- 
cantly higher stresses and larger displacements of the plate at short 
times after impact than similar sphere and cylinder projectiles. 
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XII. APPENDIX 


EVALUATION OF INTEGRALS 

In this thesis the Bateman Manuscript Project Volumes, Tables of 
Integral Transforms (ref. 6), enabled the evaluation of many difficult 
integrals that probably would have been intractable had such a set of 
volumes not been available to the author. 


Integral Occurring in the General Stress 
Expression (eq. (24)) 

The double-integral expression for the stress at the origin was 
found to be 



pIo(pu)si n (p 2T )dp 


(Al) 


The infinite integral can be evaluated from Weber's first exponential 
integral (ref. 7> P» 39^ ) 


poo 

J pj 0 ( u p) ex p(-^ 2 p 2 ) d p = 


1 ( u 2 

ZS exp [~ I^J 


2q' 


(A2) 


By letting q 2 =• It and recalling that 


exp(-iG) = cos 0 - i sin 0 


(A3) 


the integral becomes 


J ^o( u p) ^cos(p 2t ) - i sin(p2 T )jp dp = ~jsin^|^ - 


i cos 



(A4) 
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Equating real and imaginary parts yields 

J pJ 0 ( u P) c0S (p 2T ) d P ® sin ( ( A 5) 

and 

J pJ 0 (up)sin(p2 T )dp = ~ cos (a6) 

Sneddon applied this technique to Weber’s second exponential inte- 
gral and arrived at equation (ll4) on page 157 (ref. l), which reduces 
to the first Integral (eq. (A5)) by letting r go to zero. Substituting 
the second Integral (eq. (a 6)) into equation (Al) yields equation (25), 
that is, 


a-r 


• 1 h S u %( u ) cos (^ <iu 


(A7) 


Determination of gp(p) (eq. (ll)) for the 
Cone and Sphere 

Cone.- The integral 


J u 2 J 0 (pu)du 


(A8) 


which occurs in the displacement of the plate impacted by a conical 
projectile (eq. (28)) can be evaluated by changing the integral to a 
Mellin transform of 

f(u) = J 0 (pu); 0 < u < % 


= 0 


1 < u < 00 


(A9> 
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By definition the Mellin transform of a function is 


g(s) ss J f(u)u s "^ du 


(A10) 


On page 526 of reference 6, volume 1, the Mellin transform of 


f(u) » J v (pu); 0 < u < 1 

= 0 ; 1 < u < w 


(All) 


is given as 


V (l±* } v + + 1 ,- i A 

2 v (s + v)r(v + 1) ^ ^ \ ^ ^ ^ / 


(A12) 


which yields the desired integral if v is zero and s is 3# that is, 




(A13) 


Sphere . - The displacement c ' the plate impacted by a sphere con- 
tains the integral (eq. (30)) 


g sph (p) = J U ( 1 ■ u2 ) 1//£ J Q (pu)du 


(Al4) 


which can also be evaluated as a Mellin transform. On page 327 of 
reference 6, volume 1, the Mellin transform of 
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f(u) «(l - u 2 )* J v (pu); 0 < u < 1 

» 0 ; 1 < u < a? 


(A15) 


is given as 


P V B(A 


I (a 


Vi T 8 7 ^ + 1,S ~^ + 1 + x,m * p2 ) 

2 V+1 r(v + i) x 7 


provided 


and 


ReA > -1 


Res > -Rev 


(Al6) 


By letting v be zero, s be 2, and A be l/2, the desired Integral 
is found to be 


J uf 1 - u s ) 1/,S J 0 (pu)du = i £ p2 


(Air) 


The identical numerator and denominator parameters of the hypergeometric 
function cancel, leaving P 2 ^ • This hypergeometric function 


is related to the Bessel functions by the relationship 


(l p )' 


J * (p) = F(77I) ° F i( v + 1; ‘ * **) 


(A18) 


(ref. T» p. 100) 

Using this relationship, the hypergeometric function can be expressed as 

(A19) 
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Hence, the final result is 

J u(2L - u 2 ) X//fc J Q (pu)du » P“^ 2 J 5/2^ p ^ ( A2 °) 

Short Time Defl e ctions of the Plate (t « l) 

p 

For very small time after the impact, the sin p t in the inte- 
grands of the deflection integrals (eqs. (27), (29), and (51)) can be 
approximated by p 2 T to obtain the short time deflections of the plate. 

Cylinder .- The approximation in the cylinder case gives the short 
time deflection of the plate as 

n oo 

« !) * T J J^pJj^pTjJdp 

= Tt)" 1 / 2 f p' 1 / 2 J 1 (p)J 0 (pri)(ptl) 1 / 2 dp (A 21 ) 
J o 

By definition this integral is the zeroth order Hankel transform of the 
function p-l/2 Jj^p). The transform is given on page l4 of reference 6, 
volume 2, as 


f 00 -l/2 1/2 f W 2 J 0 < r\ < 1 

/ p 1/2 J 1 (p)J 0 (pn)(pn ) 1/2 dp * l 

J o ; 1 < r < 00 


(A22) 


The short time deflection of the plate for the cylindrical projectile is 


w cyl (n> T « 1) = 


TJ 0 < T) < 1 

0; 1 < T) < 00 


(A23) 
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Cone .- For the conical projectile, the short time deflection 
expression is 


00 

Ji(p)J 0 (*n)ap 


1 

-j T 


P l F 2(i ;1 4 J " F p2 ) J 0 ( pT ^ dp (A24) 


The first integral is the same as in the cylinder case; the second inte- 
gral can be rewritten as 


n_1 ^ 2 1 “ i + p 2 ) J o ^ pTl ^ pT1 ^ :1 ' //2 dp ( A 25 ) 


which by definition is the zeroth order Hankel transform of the function 
d 1 / 2 .rJiil.il . i t>2 


i F 2\2^’a’ 


p 2 ). On page 88 of reference 6, volume 2, the vth 


order Hankel transform of 



is given as 


r(y + i)r(p) g y+i 2a - v -| (1 
r(a)r(p - a) 


p.. P-a-1 

t r) ; 


0 < TJ < 1 


0 ; 1 < T) < 00 


(A26) 


(A27) 


provided Rev > -1 and Rep > Rea > 0. 
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Since the order in which the numerator and denominator parameters 
of a hypergeometric function are written makes no difference, the choice 
of v » 0, a a ”, and p ® \ in the above transform yields the desired 

<5 C 

integral, that is, 





0 < TJ < 1 

1 < TJ < oo 

(A28) 


This Integral plus the Integral evaluated for the cylinder case combine 
(eq. (A24)) to give the short time deflection for a plate impacted by a 
conical projectile as 


* co (n>T « 1) 


f t(i - n); o < n < i 

(0 } 1 < Tj < oo 


(AS?) 


Sphere . - Using the short time approximation as before, the deflec- 
tion for a plate impacted by a sphere is 


^sph( T l' T << 



J 5 / 2 ( p ) J o( pT l) dp 


(A50) 


Again rewriting this as a zeroth order Hankel transform 


”sph( T '> T << ^ = (l) 1 ^ 2 TT \" 1 ^ 2 f P" 1 J 3 / 2 ( p )J 0 ( p n)( p T)) 1 ^ 2 dp 


the integral can be evaluated from the vth order Hankel transform of 


(A51) 


the function 
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given as 


v 

P 



g v-n+i n y+§ 

r(|i - v)a^ 




0 < ti < a 


0; a < Tj < oo 


(A32) 


(A35) 


on page kQ of reference 6, volume 2, provided a > 0 and -1 < Rev < Refi. 

3 

In this case v « 0, |i * and a * 1. The desired Integral is 

p* > f(|) tj3yS(i - t j 2 ) 1 ^ 2 ; 0 < tj < 1 

J o 0 0 ;l<n<“ 


(AJ4) 


which gives 


«, , , , t(1 - T) 2 ) 1 ^ 2 ; 0 < T) < X 

v S ph(n>T « i) = 

1 0 ; 1 < i) < ■» 


(A55) 


for the short time deflection of a plate after being impacted by a 
spherical projectile. 

Maximum Stress in the Sphere Case (eq. (58a)) 

The Integral for the stress at the axis of symmetry for the sphere 

case is 

“eph ‘ & f Q u(l ' u2)l/2 ooa (^) du ( A ? 6 ) 
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Again; this integral can be handled as a Mellin transform if a change of 

o 1 

variables is made. Let x * u c and b a ^ , then 

®sph * & f ^ " x ^ 2 cos vbx)dx (A57) 

On page 320 of reference 6, volume 1, the Mellin transform of 


v»l 

f(x) a (1 - x) cos(bx)j 0 < X < 1 
a 0 \ 1 < X < oo 

Rev > 0 


(AJ8) 


is given as 


- B(s,v)j^F 1 (sjS + v;ib) + jFjjsjs + v;-lb)J 


(AJ9) 


for Res > 0. By letting v be ^ and s be 1, the desired integral 


becomes 


f o (1 - x)1/2 coB ( bx)d * - j[i r i( i >l> lb ) + i F i( i 4 i ' lb ) 


The confluent hypergeometric function is the power series 


(a4o) 


l F ^ a ’ b ’ Z ) “ 1 + b IT + b[b + l5 2! + ' ‘ ' 


V < a >n z a 
’ n =0 (b)n ” ! 


(A4l) 


4o 


where Pochammer's notation used for the numerator t md denominate* 
parameters Is 


and 


(a) n * a(a + l)(a + 2) . . . (a + n - 1) 


(a) Q * 1 


(A42) 


Adding the confluent hypergeometrlc functions of equation (A2J5) term by 
term eliminates all the imaginary terms. In general, then 


^(ajbjlz) + ^(ajbj-lz) = 2 ^ (-1)“ 


2n 


n*0 


(b) 2n (2n)! 


These results give 


(A43) 


fV 

cr 


sph 


I- Y (_l) n ^2n iitlll 

“A \1>" 


(AW) 


which reduces to ( eq. (38b)) 


suh - i n l/£ T LA. 

sph k* L ( a 

n ”0 + |j 


(A45) 
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